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ABSTRACT 

We study the relative importance of various 
damping parameters in ultraviolet lines under con- 
ditions of normal stellar atmospheres 
(10 , 000° K <_ T f £ 30 , 000°K, log g = 4) . It is 

found that the damping by electron collisions is 
always much smaller than either the radiative or 
the classical damping, except for lines that arise 
from high-excitation levels and do not therefore 
develop strong damping wings. A damping constant 
equal to 10 times the classical value is always 
much too large when compared to the more realistic 
(electron + radiative) damping. 


I . INTRODUCTION 

The spectra of 0 and B stars show a large number of metal- 
lic lines in the ultraviolet region of the spectrum — that 
is, between the Lyman limit and approximately 2000 A. We 
refer the reader to the results presented at the IAU Symposium 
36 (Houziaux and Butler, 1970) . Some of these lines are so 
strong that they can be observed with a resolution as low as 
10 A (Code and Bless, 1970). These strong lines can be heavi- 
ly saturated and hence develop important damping wings. So 
one must pay some attention to the damping parameters that 
one is going to use in line-profile calculations. In addition 
to the ordinary spectral analysis of each individual line, 
calculation of the wings of these strong lines must be carried 
out carefully because (a) they may account for a nonnegligible 
fraction of the total blanketing effect, and (b) they may 
block out a substantial amount of energy from the spectral 
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ranges (narrow or broad) to which they belong. Point (a) con- 
cerns the structure of atmospheric models, and (b) relates to 
the effects of these lines on such observable quantities as 
low-resolution spectroscopy, or narrow- to-broad-band photome- 
try. Earlier attempts to calculate blanketed models with 
ultraviolet lines have included a damping constant equal to 
10 times its classical value (see Adams and Morton, 1968; 
Hickok and Morton, 1968; Bradley and Morton, 1969; Wayne van 
Citters and Morton, 1970) . The reason is perhaps that their 
value represents empirically rather well the damping in solar- 
type atmospheres. But under solar conditions, the broadening 
by collisions with neutral particles is a predominant damping 
mechanism, whereas in the hotter atmospheres studied here, it 
is not. In this paper, we compare various damping parameters 
and from them calculate a number of line profiles. The range 
of temperatures considered is such that damping by collisions 
with neutral hydrogen can be neglected — that is, the temper- 
atures are representative of model atmospheres with effective 
temperature T e ff larger than 10,000°K. The upper limit of the 
temperature range has been taken equal to T e ff = 30,000°K 
because in hotter atmospheres, the effects of non-LTE are ex- 
pected to be very large (Mihalas and Auer, 1970) . 


II. DAMPING PARAMETERS 

For metallic lines, which are the only ones considered in 
this study, we shall consider two damping mechanisms: 

1. Natural, or radiative , damping y R , which, for a line 
with upper level u and lower level £ , is given by 


Yr — 2 A Jti + 2] ^uj * (D 

E i <E £ E j <E u 

Here, the Aj^j_'s are the transition probabilities from level £ 
to level i, and A u j , from level u to level j, measured in 
sec -1 . The summations are carried over all levels i lower 
than level £, and all levels j lower than u. The line-opacity 
profile is Lorentzian, with damping constant yr. 

2. Broadening by electron collisions. According to Sahal- 
Brdchot (1969a, b) and to Chapelle and Sahal-Brechot (1970), 
the impact approximation is valid for electron collisions 
under the conditions of normal stellar atmospheres. Hence, 
the line-opacity profile is Lorentzian, with a damping parame- 
ter Y e (sec -1 ) , where Ye is proportional to the number density 
of the perturbers, that is to the electron density Ne (cm -3 ). 



BROADENING PARAMETERS FOR ULTRAVIOLET LINE PROFILES 


379 


Because both these broadening mechanisms give rise to Lo- 
rentzian profiles, the resulting profile is Lorentzian, with 
damping y equal to 


Y = YR + Y e (sec M • ( 2 ) 

The damping due to collisions with protons or He + ions is 
much less than the electron damping (Sahal-Brechot and Segre , 
1970). As we shall see, the electron damping itself plays a 
minor role in all cases where damping wings are strong. 

Hence, we neglect the damping due to collisions with H+ and 
He + ions . We shall also compare the radiative and electron 
damping to the classical damping y c , 

Yc = 0.221 x 10 1 6 /A 2 (sec -1 ) , (3) 

where X is the wavelength of the line in angstroms. A number 
of model atmospheres with ultraviolet line blanketing have 
been computed with 10 times y c , the value of Morton and his 
coworkers cited earlier. Because the classical value for 
damping is easy to compute, it would be useful to derive an 
empirical value of the damping constant as a multiple of y Cl 
for the following reasons. First, the computation of line- 
blanketing effects involves lengthy calculations (hundreds or 
thousands of lines) ; it is desirable to keep the number within 
reasonable limits. Second, for many of these lines, the data 
for Yr and y e are simply unavailable. If, however, the radi- 
ative damping y R and the electron damping y e are available for 
the very strong lines, we can make the following predictions. 

If y R is predominant, the damping mechanism does not give any 
information on the properties of the atmosphere, because y R 
is essentially an intrinsic atomic parameter. On the other 
hand, if y e is predominant, the damping part of the line pro- 
files would give us information on the state of the gas, 
mainly the electron density and only slightly, the temperature. 
This information would then play the role of a gravity 
indicator . 

To carry out these various comparisons, we selected those 
multiplets for which the necessary data were available. These 
data, together with some other basic line data, are displayed 
in Table 1. For each multiplet, we indicate the name of the 
ion in the first column, followed by the abundance by number 
log e (with log e(H) = 12), the ionization potential X(in ev) , 
the wavelengths X (A) , and the oscillator strengths log gf of 
each multiplet' s lines. Next, we give the excitation poten- 
tials (in ev) of the lower (X^) and upper (X u ) levels. Final- 
ly, the last three columns contain the electron damping 
C e = 10 6 y e /N e at T = 20,000°K, the radiative damping yr, and 


20 ,000°K 
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the classical damping Y c . The two latter are in units 

10“ 8 sec -1 , and Y c has been calculated according to equation 

(3) . 

The values of the electron damping have been taken from 
Sahal-Brdchot and Segre (1970) . It should be noted that the 
temperature dependence of Y e is weak. In most cases of inter- 
est, Y e varies like T -0 * 25 . Also, the damping Y e is the same 
for all lines of a given multiplet (Chapelle and Sahal- 
Brechot, 1970) . 

The radiative damping Yr has been calculated according to 
equation (1) by searching through the tables of atomic-transi- 
tion probabilities by Wiese, Smith and Miles (1969) and by 
Wiese, Smith and Glennon (1966). The log gf values have also 
been taken from these tables. No critical discussion of the 
contents of these tables has been attempted. It could be 
that lifetimes of some levels are missing, which would lead 
to radiative-damping constants that are too small. (This may 
be the case of multiplet Si II at A = 1808 A.) We give only 
one value of Yr for each multiplet because the variations of 
Yr within a given multiplet are at most ±0.1 * 10 8 (sec -1 ). 

Before we discuss these various damping parameters, we 
describe Figure 1. It shows the damping log Y e and log Yr for 
each multiplet of Table 1 (in units sec -1 ). The open circles 
indicate the electron damping, and the crosses the radiative 
damping. Each multiplet is labeled by its wavelength. The 
Ye is for an electron density N e = 10 14 (cm -3 ) and a tempera- 
ture T = 20,000°K. This density corresponds to rather high 
values of the optical depth ( t 5q 00 ~ 1) i n atmospheres with 
a gravity log g = 4 and T e ff > 10,000°K. In shallower layers, 
where the lines are more likely to be formed, the electron 
damping would be lower. And in atmospheres with lower gravi- 
ties, it would be still lower. Thus, the values of Y e in Fig- 
ure 1 must be considered as upper limits in most cases of 
astrophysical interest. Also shown in Figure 1 is the classi- 
cal damping Y c , together with 10 times Y c , as a function of 
wavelength . 

From Table 1 and Figure 1, we can easily draw the follow- 
ing conclusions: 

1. The electron damping is nearly always lower than the 
radiative damping, even with the rather high value of 

N e = 10 14 (cm -3 ). The only lines for which the electron damping 
is important are those arising between very highly excited 
levels. Hence, these lines are not strong, and damping does 
not play an important role. 

2. For most lines whose lower level is not too highly ex- 
cited (and which may therefore be very strong) , the radiative 
damping is between 1/3 and 3 times the classical damping. One 
exception is the already mentioned multiplet of Si II at 
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Figure 1. — Vamping log y o ^ a numbea o ultraviolet multiples. 
Cro-ixeA Indicate radiative damping. Open circle <6 Indicate 
electron damping. The data oi each multlplet are labeled 
with lti> wavelength. Solid line* t>how classical and 10 
time -6 claxxlcal damping. 

X = 1808 A. 

3. For all lines that qualify to develop strong damping 
wings, the value Y = 10 Y c is much too large, by a factor 3 at 
least, with respect to either YR or y = Yr + Y e . 

4. The radiative damping being predominant, we do not ex- 
pect any gravity effect on the line profiles via the damping 
mechanism. This answers the question that we raised previous- 

iy. 


III. ULTRAVIOLET LINE PROFILES AND EQUIVALENT WIDTHS 

In this section, we consider how our various damping para- 
meters influence the profiles and equivalent widths of strong 
lines . These lines are supposed to be formed in stellar 
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atmospheres with effective temperatures T e ff ranging from 
10,000 to 30,000°K and with a surface gravity log g = 4 . How- 
ever, because we want to study the relative effects of the 
damping parameters with respect to each other , we do not worry 
too much about the degree of realism of our computations (i.e., 
abundances, microturbulence, etc.). So we simply indicate how 
the calculations have been performed, without any critical 
discussion of matters not directly related to the damping 
problem. 

To calculate our line profiles, we used the program ATLAS 
(Kurucz, 1969), modified to calculate metallic-line profiles 
and equivalent widths (Peytremann, 1970) . The basic assump- 
tions are those of classical model atmospheres: plane- 

parallel and in hydrostatic and radiative equilibrium. Occu- 
pation numbers are given by the Saha-Boltzmann relations. 

The source function is calculated according to the Milne- 
Eddington relation, with scattering by free electrons. (Ray- 
leigh scattering is negligible at these high temperatures.) 

The line opacities are treated in pure absorption. With these 
assumptions, it is clear that the line cores are incorrect, 
but they are rather narrow as compared to the strong damping 
wings. As stated previously, we are currently interested in 
these damping wings, and it is precisely these that contribute 
mostly to the equivalent widths of strong lines, as we shall 
see later. 

The line-opacity profiles are given by Voigt functions cal- 
culated according to Baschek, Holweger, and Traving (1966). 

At each frequency, we sum up the continuous opacities and the 
opacities of all the relevant lines (blends) . In our case, 
we take into account all lines of a given multiplet. Also, 
for clarity, we do not include the opacities of nearby hydro- 
gen lines . The equivalent widths have been computed accord- 
ing to Baschek et al. (1966) with the following exception. 

Most of the lines studied here are so broad that the residual 
fluxes 


r A = F(A)/F C (A) (4) 

are calculated with respect to the continuum F C (A) and not, as 
usual, to the continuum at the line-center wavelength (F(a) is 
the flux in the line, and F c (a) is the continuum flux). How- 
ever, since the computation of the flux at each frequency 
point is time consuming, we calculate only three continuum 
points for each multiplet and obtain the F c (a) values by para- 
bolic interpolation. This completes our description of the 
computational method. 

The actual data entering the computations are (a) the model 
atmospheres and (b) the line data. (a) The models are hydro- 
gen line-blanketed models calculated with the unmodified pro- 
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gram ATLAS (Kurucz, 1969). The ratio of helium to hydrogen 
by number is equal to 1/9 . At our temperature 
(10,000 >. T e ff 30,000°K), metals are unimportant as electron 
donors. The gravity has been taken as log g = 4. As we shall 
see, the effect of the electron density (and hence of the 
gravity) is so weak that no attempt has been made to calculate 
line profiles at lower gravities, where the contribution of 
the electron broadening would be still weaker. Finally, the 
ad hoc microturbulence parameter that enters the Voigt func- 
tion has been set equal to 5 km sec -1 . The strong lines 
studied here depend very little on the microturbulence veloci- 
ty, because their damping wings are far broader than the 
Doppler core, and some calculations with velocities equal to 
0 and 10 km sec -1 have indeed confirmed this. (b) The line 
data are contained in Table 1. The actual values of the abun- 
dances and of the partition functions are those contained in 
our version of the program ATLAS (Peytremann, 1970) . 


IV. RESULTS AND DISCUSSION 

We did not calculate the profiles of all the lines listed 
in Table 1, at each temperature. In Figures 2 to 9, we dis- 
play typical results. At the present time, we cannot compare 
our results to observations, because, for one reason, the 
instrumental profiles must be known in order to fold them in- 
to the theoretical profiles. Moreover, the definition of ob- 
served residual fluxes or equivalent widths depends on an 
accurate definition of the continuum. This is difficult, and 
perhaps impossible, if the observed spectral feature is 
blended with a large number of adjacent lines; this indeed is 
what is observed on ultraviolet spectral scans (Houziaux and 
Butler , 1970 ) . 

Figures 2 and 3 show the profiles of the multiplets 
11083 (N II), and A1206.51 (Si ill). The ordinates are the 
residual fluxes r^ defined by (4) , and the abscissas give the 
wavelength in angstroms. The profiles shown here are for 
models with T e ff = 20,000°K. The various profiles are labeled 
with y e for electron damping, Yr + Ye for electron plus radia- 
tive damping, and Yc for classical damping (Figure 3 only) . 
Also shown are the profiles with 10 times the classical damp- 
ing. Figure 3 also shows the pure Doppler profile (y = 0). 

It is seen that 10 Yc gives lines much too broad, as compared 
to the more realistic y = Yr + Ye case; Ye alone gives too 
weak profiles, but its effect is important if we compare it to 
the zero-damping case (y = 0 , Figure 3) . Figure 2, which dis- 
plays a six-component multiplet, shows a striking difference 
between the cases Ye an< ^ (Y e + Yr) and the case 10 Yc- Where- 
as the former profiles are resolved into four components, the 
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Figure 2. — ?fiO{ \llet> o i the multlplet oi N II at X - 1083 A, 
ion. a model with Tw/ = 2 0,000°K. y e , y R + Y e , and 10 y 
label pnoille-i calculated with electnon, electnon plu6 
n.adlatlve, and 10 tlmet> cla&blcal damping, n.eA pectlv ely . 



Flgun.e 3. — Same at> f Igune 2, ion. a SI III line at 1 206 . 5 1 A. 
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latter does not show any significant resolution of the multi- 
plet's components. To check this, one would need instrumental 
resolutions on the order of 0 . 5 A at A - 1100 A, provided that 
the multiplet is not blended with other lines. 

Figures 4 to 9 show the equivalent widths log W as a func- 
tion of the effective temperature log T e ££ of the models. The 
equivalent widths are in angstroms and are defined as usual by 


w = / (1 “ r x )dA / (5) 

o 

where r^ has been defined by (4) . The various curves are 
labeled with the same damping parameters given in Figures 2 
and 3 and defined in Section 2 . Each figure shows the results 
for a multiplet defined by its wavelength. Also shown are the 
ion and the number of components of each multiplet (for 
details, see Table 1). The multiplets of our sample all de- 
velop strong damping wings (compare with the pure Doppler 
case, Y = 0) . Some of them (A1083, Figure 4; A1108, Figure 5; 
A1206, Figure 6; A1334, Figure 8) show these strong damping 
wings over larcre effective temperature ranges. We reach the 
same conclusions as in §11 and from the analysis of Figures 2 
and 3. We note that although the electron damping has little 
effect, its relative contribution to the total damping Yr + Y e 
is larger at T e ff = 10, 000° K than at higher temperatures (see, 
for example. Figures 5 and 6) because in the layers at which 
the wings of these lines are formed, the electron densities 
are higher at T e ff = 10,000 °K, than at higher temperatures. 

An indirect consequence is that if any non-LTE effects exist 
in these regions of line-wing formation, they would be en- 
hanced at higher temperatures . 

We now summarize. As stated before, the conditions consid- 
ered here are those that prevail in atmospheres hotter than, 
say, 10,000 °K, that is, where collision broadening by neutral 
particles is negligible. 

First a damping constant of 10 times the classical value is 
much too large, by comparison with the more realistic case of 
electron-plus-radiative damping. This is not necessarily true 
for lines whose lower level is highly excited, but these lines 
are not expected to form strong damping wings anyway. 

Second, for the strong lines we are dealing with, the ef- 
fect of electron damping is generally much smaller than that 
of radiative damping. So we do not expect observable gravity 
effects due to the broadening mechanisms. 

Third, in most cases studied here, the electron-plus-radia- 
tive damping gives more or less the same results as the clas- 
sical damping. Compared to many other sources of errors, 
these differences are not dramatic. For calculations that 
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involve hundreds or thousands of lines (blanketed model atmos- 
pheres) or that concern lines with unavailable damping data, 
we would therefore suggest, as a rule of thumb, to use the 
classical damping constant rather than that 10 times greater . 

This work has been done while the author held an ESRO-NASA 
International Fellowship. The research was supported in part 
by contract NAS 5-1535 from the National Aeronautics and 
Space Administration. 
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